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In this paper we systematically study a simple class of translation-symmetry protected topological 
orders in quantum spin systems using slave-particle approach. The spin systems on square lattice 
are translation invariant, but may break any other symmetries. We consider topologically ordered 
ground states that do not spontaneously break any symmetry. Those states can be described by Z2A 
or Z2B projective symmetry group. We find that the Z2A translation symmetric topological orders 
can still be divided into 16 sub-classes corresponding to 16 new translation-symmetry protected 
topological orders. We introduced four Z2 topological indices = 0, 1 at k = (0,0), (0,7r) , (n,0), 
(n, tt) to characterize those 16 new topological orders. We calculated the topological degeneracies 
and crystal momenta for those 16 topological phases on even-by-even, even-by-odd, odd-by-even, 
and odd-by-odd lattices, which allows us to physically measure such topological orders. We predict 
the appearance of gapless fermionic excitations at the quantum phase transitions between those 
symmetry protected topological orders. Our result can be generalized to any dimensions. We find 
256 translation-symmetry protected Z2A topological orders for a system on 3D lattice. 
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I. INTRODUCTION AND MOTIVATION 

For long time, people believe that Landau symme- 
try breaking theory 1 and the associated local order 
parameters 2 ' 3 describe all kinds of phases and phase 
transitions. However, in last 20 years, it became more 
and more clear that Landau theory cannot describe all 
quantum states of matter (the states of matter at zero 
temperature). 4,5 A nontrivial example of new states of 
matter beyond Laudau theory is the fractional quantum 
hall (FQH) states. 6 FQH states do not break any sym- 
metry and hence cannot be described by broken sym- 
metries. The subtle structures that distinguish different 
FQH states are called topological order. 4,5,7 Physically, 
topological order describes the internal order (or more 
precisely, the long range entanglement) in a gapped quan- 
tum ground state. It can be (partially) characterized by 
robust ground-state degeneracy. 4 ' 8 Recently, many dif- 
ferent systems with topologically ordered ground states 
were found. 8-16 

Quantum spin liquid states in general contain non- 
trivial topological orders. In the projective con- 
struction of spin liquid (also known as slave-particle 
approach), 17-23 there exist many spin liquids with low 
energy SU(2), U(l) or Z 2 gauge structures. Those spin 
liquids all have the exactly the same symmetry. To distin- 
guish those spin liquids, we note that although the spin 
liquids have the same symmetry, within the projective 
construction, their ansatz are not directly invariant un- 
der the translations. The ansatz are invariant under the 
translations followed by different gauge transformations. 
So the invariant group of the ansatz are different. We 
can use the invariant group of the ansatz to characterize 



the new order in the spin liquids. The invariant group of 
an ansatz is formed by all the combined symmetry trans- 
formations and the gauge transformations that leave the 
ansatz invariant. Such a group is called the Projective 
Symmetry Group (PSG). 24 Thus although one cannot 
use symmetry and order parameter to describe different 
orders in the spin liquids, one can use the PSG to charac- 
terize/distinguish the different quantum/topological or- 
ders of spin liquid states. 

The simplest kind of topological orders is the Z 2 topo- 
logical order where the slave-particle ansatz is invariant 
under a Z 2 gauge transformation. According to the PSG 
characterization within the projective construction, for 
system with only lattice translation symmetry, there can 
be two different classes of Z 2 topological orders labeled 
by Z2A and Z2B (a Z2B ansatz has tt flux going through 
each plaquette). 24 In this paper, we will study the Z2A 
topological orders and ask "are there distinct Z2A topo- 
logical orders?" We find that there are indeed distinct 
Z2A topological orders. They can be labeled by four 
Z 2 topological indices Cfc = 0, 1 at k = (0,0), (0, tt) , 
(tt, 0), (7r,7r). So the £fc characterization is beyond the 
PSG characterization of quantum/topological order and 
provides additional information for translation-symmetry 
protected Z 2 topological order. 



II. A GENERAL "MEAN-FIELD" FERMION 
HAMILTONIAN OF Z 2 TOPOLOGICAL ORDERS 

We will use the projective construction (the slave- 
particle theory) 17,23 to systematically construct different 
translation symmetric Z 2 topological orders in spin-1/2 
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systems on square lattice. In such a construction, we 
start with "mean-field" fermion Hamiltonian 24 

#mcan = E V>j + ^(^Ivij^Pj + h.C.) + E ^i a i^i 
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(1) 

where Uij, r}ij, di are 2 by 2 complex matrices. The r\- 
term is included since our spin- 1/2 systems in general do 
not have any spin rotation symmetry. We like to mention 
that <Zj are not free parameters, ai should be chosen such 
that 

(*^W.o*)|^V^| $ (^W.a*) ) = 0j Z ^ 1,2,3 (2) 

where a 1 are the Pauli matrices. In this paper, we will 
only consider translation invariant ansatz Uij = Ui+ a ,j+a 
and rjij = rji +a _j +a . Those states are characterized by 
Z2A PSG and are Z2A topological states. 24 

Let |^ / mcan ? * J ' ) } be the ground state of H mcan . Then 
a many-spin state can be obtained from the mean-field 
state \^eln ij) ) by projection 



I* 



{Uij,T}ij)\ 

spin / 



'p|\Tr( u *rf» 7 fti)\ 



(3) 



into the subspace with even numbers of fermion per site. 
Here the projection operator is 



l + (-l) r 



and Ui = V'iV'i is fermion operator at site i. 

We note that after the projection, each site can have 
either no fermion or two fermions. If we associate the 
no-fermion state as the spin-down state and the two- 
fermion state as the spin-up state, then the projected 

state ^^l^ 11 ^) can be viewed as a quantum state for 
the spin system. This is how we construct many-spin 
state from the mean-field Hamiltonian. For each choice of 
the ansatz {uij, rjij, ai), this procedure produces a phys- 
ical many-spin wave function \ f { s l;i' Vij) ). So the ansatz 
) can also be viewed as a set of labels that label 
a many-spin state and ^phT^) can be view as a trial 
wave function for a spin-1/2 system, with and rjij 
being variational parameters. 

For the translation invariant ansatz, one can rewrite 
the "mean-field" fermion Hamiltonian in momentum 
space by introducing 



V 4,- k J 



and 



Note that satisfy the following algebra 

{*/fc. *Jk'} = $ijtk-k>, {*/*, *./fc'} = FijSk+k 
where 

Vi 



r = CTi <g> (To 



o-i 



(4) 



We also note that (*t fe , *-fc) can be expressed in term 
of 



(5) 



In terms of ^fc, H mcan can be written as 



fe/O 



fe=0 



where ~ir < k x ,k y < +ir and M(k) arc 4x4 Hermi- 
tian matrices M(k) = M^(fe). Here k = means that 
(k x ,k y ) = (0,0), (0, 7r), (7r,0), or (ir, n). Also fca; and k y 
are quantized: fc^ = j~-x integer and k y = j^-x inte- 
ger, where L x and L y are size of the square lattice in 
the x- and y-directions. Note that on an even by even 
lattice (ie L x = even and L y = even), (k x ,k y ) = (0,0), 
(0, 7r), (tt, 0), or (it, it) all satisfy the quantization condi- 
tions k x = jf-x integer and k y = jf-x integer. In this 

case, J2k=o sums over an tne f° ur points (k x , k y ) = (0, 0), 
(0,7r), (it, 0), and (tt,tt). On other lattices, ^2 k=0 sums 
over less points. Say on an odd by odd lattice, only 
(k x ,k y ) — (0,0) satisfies the quantization conditions 



k x = ^f- x integer and k y = j^-x integer. In this case, 

Sfe=o sums ° ver ° n ly (k x ,k y ) = (0,0) point. 
We note that 

*j> fc = 2, <t ® o-3*fc = 0. 

Thus up to a constant in H mcan , we may assume M(fe) 
to satisfy 

TrM(fc) = 0, Tr[M(k)a <g> cr 3 ] = 0. (7) 
Due to Eq. (5), 

*L fc M(-fc)*_ fc = TrAf(fc) - tf£rM T (-fc)r* fc . 
Thus, we may rewrite Eq. (6) as 

fe>0 fe=0 

U(k) = M(k) - TM T {-k)T. (8) 



Here k > means that fc ^ and fc y > or k y = 0, k x > 
0. Clearly U(k) satisfy 

c/(fe) = -r?/ T (-fc)r, U(k) = U\k). 

Now we expand U(k) by 16 Hermitian matrices 

M {aj3 } ee cr Q (gicr/3, a, /3 = 0,1,2,3, (9) 
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where ao = 1. We have 

U(k) = J2 c {*P}( k )M {af3} 
{«,/?} 

where C{ Q/ g} (k) are real. The 16 4x4 matrices M{ Q( g} can 
be divided into two classes : in one class, the matrices 
satisfy 

M = —TM T Y. 

We call them " even matrices"; In the other class, the 
matrices satisfy 

M = TM T Y. 

We call them " odd matrices" . 
There are six even matrices, 

-^{30} = °~3 ® CT , M{- 12 } = £7i ® £7 2 , M {22 } = <7 2 ® CT 2 , 
^{33} =CT3® CT 3 , M {31} = CT 3 ® £7i, M {02} = <7 ® <7 2 . 

For above six matrices, M{ 30 }, M{i 2 } and M{ 2 2} anti- 
commute with each other, 



{M {30} ,M {12} } = 0, 
{M {30} ,M {22} } = 0, 
{M {12} ,M {22} } = 0. 



(10) 



III. CLASSIFICATION OF Z2A TOPOLOGICAL 
ORDERS 



For a generic choice of Uij and r/ij, the corresponding 
mean-field Hamiltonian (8) is gapped. Note that the en- 
ergy levels of the mean-field Hamiltonian (8) appear in 
(E, —E) pairs. The mean field state is obtained by filling 
all the negative energy levels. The mean-field Hamil- 
tonian is gapped if the minimal positive energy level is 
finite. The gapped mean-field Hamiltonian corresponds 
to a gapped Z2A spin liquid. 

As we change the mean- field parameters Uij and rjij, 
the mean-field energy gap (and the correspond energy 
gap for the Z2A spin liquid) may close which indicates 
a quantum phase transition. Thus if two gapped regions 
are always separated by a gapless region, then the two 
gapped regions will correspond to two different quantum 
phases. We may say that the two quantum phases carry 
different topological orders. 

In the following, we introduce topological indices that 
can be calculated for each gapped mean-field ansatz 
{ u ij,Vij)- We will show that two gapped mean- field 
ansatz with different topological indices cannot smoothly 
deform into each other without closing the energy gap. 
Therefore, the topological indices characterize different 
Z2A topological orders with translation symmetry. 



M{ 33 }, M{ 31 jand M{ 2} anti-commute each other 



{M {33} ,M {31} } = 0, 
{M {33} ,M {02} } = 0, 
{M {31} ,M {02} } = 0. 



While each of M. 



M(, 2) and M. 



{30} 7 Jw {12} 



{22} 



(11) 

commutate with 

each of M{ 33 j,M{ 31 } and M{ 2}- For the coefficients of 
the even matrices C{ Q( 3}(fc) in the "mean-field" fermion 
Hamiltonian, {a[3} = {30}, {12}, {22}, {33}, {31}, {02}, 
we have 

C{a/3}(fc) = C{ Q(3 }(-fc)- 

In addition, there exist ten odd matrices: 



M. 



{00} — ^0 
-^{03} = CO 



M. 



{23} — a 2 
M{01} = (To 



' 0-3, 
' OS, 

' 0-1. 



M{ W } = (71 <g> <7 , 
M { n } = (71 ® £71, 
-^{32} = O-3 <g> £7 2 , 



M. 
M 
Mi 



{20} 
{13} 



(7 2 
(71 



{21} — °~2 



1 c , 

1 03, 
1 Cl, 



For the coefficients of odd matrices C{ Q(3 j(fc) in the 
"mean-field" fermion Hamiltonian, we have 

c {af3}{k) = -C{ Q(3 }(-fc). 

Thus for odd matrices, C{ a/ 3}(fc) are odd functions of 
k Xl k y and are fixed to be zero at momentum (0, 0), (0, n), 
(tt,0), (7T,7r) 



A. The topological indices 

In the section II, we obtained the "mean-field" fermion 
Hamiltonian in momentum space (8), which has a form 
H mcan = H(k > 0) moan + H(k = 0) moan . Let us di- 
agonalizing the "mean-field" fermion Hamiltonian at the 
points k > as an example. Introducing 



W(k)V k = 



( a k \ 



(12) 



where 



W{k)U(k)W\k) 

ei(fe) 

-£i(fe) 

£ 2 (fe) 

-£ 2 (fe) 



(13) 



C{ a p}(k = 0) = 0. 



£i(fe) > 0, and £ 2 (fc) > 0, we find 



4 



H(k > 0) moan = ^2 [ £ i( k ) a k a k + £2(k)fiif3 k - £i(fc)a_ fc aL fc - e 2 {k)(3_ k t _ k 

fc>0 

We note that both a± k and [5± k will annihilate the mean-field ground state l^mean), 

a±fe|*mcan) = 0, /3±fc|*mcan) = 0. 

At the four k — points, only the even M a ^ appear and we diagonalized the Hamiltonian differently. The four 
eigenvalues of U(k) are given by 



and 



s±±{k) = ±Jc 2 30} (k) + c\ l2} {k) + c\ 22} {k) ± Jc 2 33} (k) + c 2 31} (k) + c 2 02} (k) 



H(k = 0) mcan = i [ £ ++( k ) a k a k + £ +-{k)Pifik + s--(k)a k a f k + e_ + (fc)/? fe /?t 



k=0 



^2 [e++(k)ala k + e+-(k)p{p k 



Const. 



(14) 



where 



W{k)U{k)W\k) 

/e++(fc) 

-e++{k) 



V 






£+-(fc) 
-£_, 



(15) 



We note that W{k) diagonalizes the linear combination 
of M {30} , M {12 }, and M {22 } in H me&n : 



W (c{ 30 }M {30 } + C{12}M{ 12 } + C{22}M{22}) VK f 



\/ C {30} + C {12} + C {22 ! ■ ' ' i : « ' ' ■ 



(16) 



W(k) also diagonalizes the linear combination of M{ 33 }, 
M{ 31 }, and M {02 } in H mcan : 

W (c {33} M {33} + c {31} M {31} + c {02 }M {02} ) IT 1 

= \/ C {33}+ C {31}+ C {02} M {33}- (17) 

So W^(fc) changes M 33 = tr 3 ® cr 3 to 

= a(fc)cr 3 ® cr 3 + b(k)a ® <ri + c(fc)a 3 ® cr 2 . (18) 

where a 2 (fc) + 6 2 (fc) + c 2 (fe) = 1. 

The energy spectrum at fc = motivates us to intro- 
duce Cfc as the four topological indices, one for each k = 
point: 



a = i-e[ e+ _(fc)]. 



(19) 



®[ C {30}(^) + C {12}( k ) + c {22}(^) 
{33} W - c {31}(k) - c {02}(^)] 



where 8(x) = 1 if x > and 6(2;) = if a; < 0. If two 
topological ordered states have different sets of topo- 
logical indices (Cfe=(o,o), Cfe=(f,o)> Cfe=(o,7r), 
then as we deform one state smoothly into 
the other, some ( k must change. ( k can only 

change when c 2 {30} (k) + c 2 {12} (k) + c 2 {22} {k) 



{33} 



(*0 



1 i(A;) + C| 02 |(fc). At that point, the 



-{3i}* 



topological ordered state becomes gapless indicating 
a quantum phase transition. Therefore, there are 16 
different translation invariant Z2A spin liquids labeled 

by Cfc=(0.0)> Cfe=(0,7r); Cfc=(7r.0)j Cfc=(7r,7r) = HH, 1100, 

1010, 1001, 0101, 0011, 0110, 0000, 1000, 0100, 0010, 
0001, 1110, 1101, 1011, 0111. 



B. The topological degeneracy 

Let's calculate topological degeneracies for different 
topological orders through the projected construction. 



(u\ 



Now we use \m,n) = |*moin "' iS ) (|m,n) = |0,0), 
|0, 1), 1 1, 0) , |1,1)) to note four degenerate ground 

states on a torus. Here (u^' n \ ^™'"^) is defined 
as ((_)™^i)(_)™w(y) Uii) (-)™*(ij)(-)™y(ij) mj ). 
s x ,y(ij) have values or 1, with s XtV (ij) = 1 if the 
link ij crosses the x — L x , or y = L y line(s) and 
s x, y (ij) = otherwise 24 . In fact, the four mean-field 

states |m,n) = l^melii 13 ) are obtained by giving 
the fermion wave- functions ip(x,y) different boundary 
conditions: 



— c 



1>(x,y) = (-l) m i/>(x, y + L y ), 
ip(x,y) = {-l) n il>(x + L x , y). 



(20) 



■5 



To obtain a physical state from the mean-field ansatz 

l^mein^}, one needs to project \^Hn ij) ) into the sub- 
space with even numbers of tp fermion per site. So the 
mean-field state must have even numbers of tp fermions 
in order for the projection to be non-zero. The total tp 
fermion number has a form N = N k ^ + ^fe=o where 

N k ^o = E (V'feV'fe + i>- k ip-k) and N k=0 = J2 V'feV'fe- 

fc>0 fc=0 

We note that, for k > 0, 

^_^V'i ifc V'i,fc+V'l ! _feV'i,-fc+V'2,fc 1 /'2,fc+V'| i _j s V'2,-fe 

= (_^l,fc^l.fc-V'l,-fc^l,_|,+V , I, fc ^2,fc-V , 2,-fcV'J i _ fc 



= (_)«l«fc+/3fc/3 fc +aL fc «-fc+/3^ fc /3- fc 

Hence we have 

(V"! l.V'a.k+^l -fcV'a.-fc) 



(21) 



(_).=. 
(-) 



vi>(««.'7«)\\ 
* moan // 



_ I , T/ ( U «.»)i3)\ 

— 1 mean / 



(22) 



for k > 0. We see that the total number of the tp fermions 
on all the fe ^ orbitals is always even. 

So to determine if the mean-field ground state contain 
even or odd number of tp fermions, we only need to count 
the number tp fermions at the four special points: (0, 0), 
(0,tt), (tt,0), (7r,7r). For k = 0, 

= (_)3*i^3«>ff3*fc i (23) 

Using eqn. (18), we find 

^l<7 3 ® cr 3 *fe 
= ^a(fc)(c4.a fe - a fe a|. - + 0kP k ) 
+ \Kk){a{Pk - "k i,, + P{a k - 0hOt{) 

= a(k)(a{a k - + fc^X"*/ 3 * + ^ k a k ) 

+ ic(k)(- a {l3 k + fi{ak). (24) 

We see that 

(_)^4,t>i.*+v4,j,^2.fc (25) 

_ e 7ri(a(fe)(ata fc -/3t/3 fe )+6(fe)(at/3 fc +/3t afc ) +ic ( fe) (_ a t^ fe+/3 t afc )) 

If we treat (a k ,f3 k ) as an iso-spin-1/2 doublet, then the 
above operator generates a 2tt rotation since a 2 (A;) + 



b 2 (k) + c 2 (fc) = 1. This is consistent with the follow- 
ing relation 

afe (_)V'I, fc '/'l,fe+V4 ife V>2,fc = _(_)^ t ,fc'Al,fc+'/>2, fc 'A2,fc Q , fc 

/ g fe (_)V'i t , fc '/'i,fc+V'J, fc V'2,fc _ _(_)^i t , fc </-i,fc+'/> 2 t , fc </<2,fc^ fc (26) 

Since the operator (— ) Q fc Qfc +' 3 fc' 3fc has the same algebra 
as (-^l.h^.fc+^.k^,* an d the two operators are equal 
when a(k) = 1, b(k) = c(k) = 0, we have 

(_)^!, fc , Ai,fc+'/'l, fc V'2,fc _ ^_^ot ] k a k +()lf) k _ (27) 

Thus, the total fermion number at the four points, fc = 
(0,0), (0,7r), (tt, 0), (7r, 7r) satisfies 



( _)^=o = ( _ )fc 5„ (aIafc+/3 ^ fc) . 



(28) 



For a given point of k = 0, the energies for the particles a 
and /3 are e ++ (k) and (fc). There are two situations: 

1. For e + _(fc) > 0, the fc-orbital will be filled parti- 
cle: zero a particle and zero (5 particle. 

2. For e_| (fe) < 0, the fc-orbital will be filled 1 parti- 
cle: zero a particle and one [3 particle. 

We find that the total number of tp fermion at the k = 
points and at all k are given by 



N k=0 mod 2 = N mod 2 = ^ Cfe mod 2. 



(29) 



«=o 



Note that on even- by-even lattice (ee) , all the four k = 
points k = (0,0), (0,7r), (n,0), (tt,tt) are allowed. In this 
case ./V mod 2 is the sum of all four ( k =o mod 2. On even- 
by-odd lattice (eo), only two k = points k = (0,0), 
(n, 0) arc allowed. In this case N mod 2 = C(o,o) + C(tt,o) 
mod 2. 

Let's use the topological order 1000 as an example to 
demonstrate a detailed calculation of the ground state de- 
generacy. There are four degenerate mean-field ground 

states \m,n) = |* m ean ^ ), m,n — 0,1. On 
an even-by-even (ee) lattice, the state |0,0) has peri- 
odic boundary conditions along both x and y directions. 
Among the four k = points, only k = (0, 0) point has 
Cfe = 1 as indicated by the first 1 in the label 1000. As a 
result, the ground state |0, 0} contains an odd number of 
fermions and is un-physical: P|0,0) = 0. For the states 
|0, 1), |1,0), |1,1), the k = (0,0) point is not allowed. 
Consequently, N = mod 2. Hence |0,1), |1,0), |1,1) 
are all physical states. This gave rise to three degenerate 
ground states for the topological order 1000 on an (ee) 
lattice. 

Second, we calculate the ground state degeneracy on 
an even-by-odd (eo) [or odd-by-even (oe)] lattice. For the 
state |0, 0), only two k = points are allowed: k = (0, 0) 
and k = (7r,0). Due to C(o,o) = 1) |0, 0) is forbidden, 
V\0,0) = 0. For other states |0, 1), |1,0), |1,1), k = 
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(0, 0) is not allowed on an (eo) [or (oe)] lattice. For the 
same reason, one obtains three degenerate ground states 
1 0, 1) , 1 1, 0), |1,1) on an (eo) [or (oe)] lattice. Third, we 
calculate the ground state degeneracy on an odd-by-odd 
(oo) lattice. For the state |0, 0), there is only one k = 
point: k = (0,0). As a result, |0, 0) is not permitted by 
the projection operator, P\0,0) = 0. However, without 



the point k — (0,0), other states |0, 1), |1,0), |1,1) are 
all physical. One also obtains three degenerate ground 
states 0,1), 1,0), 1,1) on an (oo) lattice. 

By this method, we obtain topological degeneracies on 
different lattices for the 16 topological orders. The results 
are given in the following table: 



1111 1110 1101 1011 0111 1100 0011 1001 0110 1010 0101 1000 0100 0010 0001 0000 

(ee) 4333344444433334 
(eo) 4333344222233334 
(oe) 4333322224433334 
(oo) .111122222233334 



C. The crystal momentum 

Next, we calculate the crystal momentum for 16 topo- 
logical ordered states. Because the spin Hamiltonian is 
translation invariancc, the ground states carry definite 
crystal momentum. To calculate the crystal momentum 
K, we note that the fermion wave function satisfies the 
(anti) periodic boundary condition. 

The crystal momentum is given by 



K 



*spin) = | *spin) 



spi 



„}• (30) 



k=0 



, (m,n) (m 



The ground state l^moin %3 ) at k ^ has a form 
ft a kCt-k(3k(3-k\0)-4> where |0)^, is the state with no ip 

k>0 

fermion. Thus, the total crystal momentum K are ob- 
tained as 



K I * s 



in) = 

k=0 



*spi„) 



(31) 



fc=0 



where we have used (at k = 0) 

V-lVfc mod 2 = Cfe. (32) 

Thus, to determine the crystal momentum, we only need 
to focus on the cases at k = 0. 

By this method, we obtain the crystal momenta on 
different lattices for 16 topological orders. The follow- 
ing tables show the crystal momenta of different ground 
states: 

K (0000) (ee) (eo) (oe) (oo) 

(00) (0,0) (0,0) (0,0) (0,0) 

(01) (0,0) (0,0) (0,0) (0,0) 

(10) (0,0) (0,0) (0,0) (0,0) 

(11) (0,0) (0,0) (0,0) (0,0) 



K (0011) 


(ee) 


(eo) 


(oe) 


(oo) 


(00) 


(0.7T) 


- 


(0,0) 


(0,0) 


(01) 


(0,0) 




(0,0) 


(0,0) 


(10) 


(0,0) 


(0,0) 


(0.7T) 




(11) 


(0,0) 


(0,0) 


(0,0) 


— 


K (1100) 


(ee) 


(eo) 


(oe) 


(oo) 


(00) 


(0.7T) 


— 


(0,7T) 


- 


(01) 


(0,0) 




(0,0) 




(10) 


(0,0) 


(0,0) 


(0,0) 


(0,0) 


(11) 


(0,0) 


(0,0) 


(0,0) 


(0,0) 


K (1111) 


(ee) 


(eo) 


(oe) 


(oo) 


(00) 


(0,0) 


(tt.O) 


(0,0 


- 


(01) 


(0,0) 


(tt.O) 


(0,0) 


- 


(10) 


(0,0) 


(0,0) 


(o,0 




(11) 


(0,0) 


(0,0) 


(0,0) 




K (0101) 


(ee) 


(eo) 


(oe) 


(oo) 


(00) 


(tt.0) 


(0,0) 




(0,0) 


(01) 


(0,0) 


(tt.O) 


(0,0) 




(10) 


(0,0) 


(0,0) 




(0,0) 


(11) 


(0,0) 


(0,0) 


(0,0) 




K (1010) 


(ee) 


(eo) 


(oe) 


(oo) 


(00) 


(tt.0) 


(tt.O) 






(01) 


(0,0) 


(0,0) 


(0,0) 


(0,0) 


(10) 


(0,0) 


(0,0) 






(11) 


(0,0) 


(0,0) 


(0,0) 


(0,0) 


K (0110) 


(ee) 


(eo) 


(oe) 


(oo) 


(00) 


(7T,7r) 






(0,0) 


(01) 


(0,0) 




(0,0) 




(10) 


(0,0) 


(0,0) 






(11) 


(0,0) 


(0,0) 


(0,0) 


(0,0) 
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K (1001) 


(ee) 


(eo) 


(oe) 


(00) 


(00) 


( 7T, 7T J 








(01) 


(0,0) 


- 


(0,0) 


(0,0) 


(10) 


(0,0) 


(0,0) 


— 


(0,0) 


(11) 


(0,0) 


(0,0) 


(0,0) 


— 


K (1000) 


(ee) 


(eo) 


(oe) 


(oo) 


(00) 










(01) 


(0,0) 


(0,0) 


(0,0) 


(0,0) 


(10) 


(0,0) 


(0,0) 


(0,0) 


(0,0) 


(11) 


(0,0) 


(0,0) 


(0,0) 


(0,0) 


K (0100) 


(ee) 


(eo) 


(oe) 


(oo) 


(00) 




(0,0) 




(0,0) 


(01) 


(0,0) 




(0,0) 




(10) 


(0,0) 


(0,0) 


(0,0) 


(0,0) 


(11) 


(0,0) 


(0,0) 


(0,0) 


(0,0) 



K (0010) 


fee) 


(eo) 


foe) 


(oo) 


(00) 






(0,0) 


(0,0) 


(01) 


(0,0) 


(0,0) 


(0,0) 


(0,0) 


(10) 


(0,0) 


(0,0) 






(11) 


(0,0) 


(0,0) 


(0,0) 


(0,0) 


K (0001) 


(ee) 


(eo) 


(oe) 


(oo) 


(00) 




(0,0) 


(0,0) 


(0,0) 


(01) 


(0,0) 




(0,0) 


(0,0) 


(10) 


(0,0) 


(0,0) 




(0,0) 


(11) 


(0,0) 


(0,0) 


(0,0) 




K (0111) 


(ee) 


(eo) 


(oe) 


(oo) 


(00) 








(0,0) 


(01) 


(0,0) 


(tt.O) 


(0,0) 




(10) 


(0,0) 


(0,0) 






(11) 


(0,0) 


(0,0) 


(0,0) 




K (1011) 


(ee) 


(eo) 


(oe) 


(oo) 


(00) 




(tt.O) 






(01) 


(0,0) 




(0,0) 


(0,0) 


(10) 


(0,0) 


(0,0) 


(0,^ 




(11) 


(0,0) 


(0,0) 


(0,0) 





K (1101) (ee) (eo) (oe) (oo) 

(00) - (0,0 - 

(01) (0,0) (tt,0) (0,0) - 

(10) (0,0) (0,0) (0,0) 

(11) (0,0) (0,0) (0,0) - 



K (1110) (ee) (eo) (oe) (oo) 

(00) - (tt,0) (0,tt) 

(01) (0,0) - (0,0) - 

(10) (0,0) (0,0) - 

(11) (0,0) (0,0) (0,0) (0,0) 



IV. HOW MANY DISTINCT Z 2 TOPOLOGICAL 
ORDERS? 

We would like to point out that the four topological 
indices at k = really describe 16 classes of mean- 
field ansatz. It is not clear if different mean-field ansatz 
give rise to different many-body spin wave functions. So 
it is possible that the 16 sets of topological indices £fc 
describe less than 16 class of Z 2 topological orders. 

On the other hand, if two Z 2 topological phases can be 
separated through measurable physical quantities, such 
as ground state degeneracy and crystal momenta, then 
the two topological phases will be really distinct. 

Using the ground state degeneracies on differ- 
ent types of lattice we can group the 16 sets 
of topological indices into 7 groups: {0000}, 
{0001,0010,0100,1000}, {0101,1010}, {1100,0011}, 
{1001,0110}, {0111,1011,1101,1110}, {1111}. Each 
group have the same ground state degeneracies. So we 
have at least 7 distinct Z 2 topological phases. If we 
assume that the boundary condition labels (to, n) are 
not physically observable, the crystal momentum distri- 
butions cannot further separate the above 7 groups into 
smaller groups. However, it is likely that the boundary 
condition labels (m, n) are physically observable by 
moving the unique type of fermionic excitations (the 
spinons) around the torus. In this case all 16 sets 
of topological indices Cfc label distinct Z 2 topological 
phases. 

V. TWO EXAMPLES OF Z2A TOPOLOGICAL 
ORDERS 

Let us discuss two different translation symmetric Z2A 
topological orders studied in 25 in more detail. The first 
one is described by the following "mean-field" fermion 
Hamiltonian within the projective construction 8 : 

■ffmcan = ^Uijlpj + ^ ^{a^i (33) 

ij i 
3 

l^i^i+x u i,i+y \0 ? 

1 i \ 2 

Ui,i+x+y = TjO + AO , 

1 \ 2 

Ui,i-x+y = TjO A<T , 

a\ = v. 

where i[> T = (ipi,^)- The other Z2A spin liquid 
comes from an exact soluble spin-1/2 model on square 
lattice - the Wen-plaquette model. 16 Its Hamiltonian 

is H = l6gJ2i S i S i+ x s i+x+y S !+y Usin g projective 
construction, one can introduce a "mean-field" fermion 
Hamiltonian: 24 to describe such a spin liquid: 

where I,J=1, 2. It is known that the ground states for 
g < and g > have different symmetry protected topo- 
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logical orders. The ground state (Z2A topological order) 
for g < is described by mean-field ansatz —rj ii+x = 
u hi+x = 1+er 3 and -r\i,i+ y = u ht+y = l-er 3 . The ground 
state (Z2B topological order) for g > is described by 
mean-field ansatz —r)i^+ x = Ui^ +X = (— ) ly (l + er 3 ) and 

'Hiji+y — Ui,i+y — 1 (J . 

We like to ask: whether the topological order for Z2A 
gapped state in Eq. (33) and the topological order for the 
Wen-plaquette model in Eq.(34) are the same one. The 
four Z 2 topological variables Cfc=o can help us to answer 
the question. 

For the Z2A gapped state, a "mean-field" fcrmion 
Hamiltonian in momentum space becomes 



where 



,{k)=Y,^ k U(k)^ k + h.c. 

k) 



(35) 



U(k) 

C{33} 
C{31} 
C{02} 



c {af3}( k ) M {a/3}, 

cos k,. 



cos k x 
r/(cos(k x 
A(cos(fc x 



+ k y ) + cos(fc x - k y )) H 
+ ky) - cos(k x - k y )). 



Because M{ 33 }, M{ 31 }.and M^ 02 y make up of an anticom- 
mutation base, we have 



e+_(fc) = - ^Jc 2 {33} (k) + c 2 31} (fc) + c 2 {02} (k) 

— yj (cos + cosky) 2 + [i](cos(k x + k y ) + cos(k x - k y )) + vf + [\(cos(k x + k y ) - cos(k x - k y ))) 2 



For the Z2A gapped state, one has 

C(0,0) = 1; C(0,tt) = 1; C(tt,0) = 1) C(7r,7r) = 1- 

It belongs to the 1111 type of the topological order. The 
topological degeneracy is 4, 4, 4, for even-by-even, even- 
by-odd, odd-by-cven respectively. On odd-by-odd lat- 
tices, the state has no energy gap. 

Another example of topological order is the Wen- 
plaquette model. The "mean-field" ansatz in momentum 
space now becomes 

U(k) = Y c {a f3}(k)M {af3} , 

C{30} = cos k x + cos ky , 

c {33} = cosk x - cosky, 

C{20} = sin fca, + sin ky, 

C{23} = sink x - sin k y . 

There are two even matrices in the ansatz, M{ 30 } and 
M| 33 }. We have 

(fc = 0) = Icosfcc + cosk v \ — \cosk x — cosfcjj . 

For the topological order of the Wen-plaquette model, 
one has 

C(0,0) = 0, C(0,tt) = 1, C(tt,0) = 1, C(7r,7r) = 0. 

It belongs to the 0110 type of the topological order. The 
topological degeneracy is 4, 2, 2, 2 for (ee), (eo), (oe), 
(oo) lattices. Because topological order in the toric-code 
model with translation invariance 12 and that in the Wen- 
plaquette model are equivalent, one can use the same 
wave-function to describe the toric code model. 



Then when one changes the Z2A gapped state denoted 
by 1111, into the topological order of the Wen-plaquette 
model denoted by 0110, quantum phase transition occurs 
with emergent massless fermion at fc = (0, 0) and k = 

(7T,7r). 

VI. NON-ABELIAN TOPOLOGICAL STATES 

We like to point out that for the four Z2A topological 
states described by {OJ = 1000, 0100, 0010, 0001, the 
corresponding mean-field Hamiltonian describes a super- 
conducting state whose band structure has an odd wind- 
ing number. 26,27 So those four Z2A topological states are 
closely related to the topological spin liquid state ob- 
tained by the projection of p x + ip y SC states. 27 As a 
result, the four Z2A topological states have a topolog- 
ical order described by Ising topological quantum field 
theory (which is the same topological quantum field the- 
ory describing the non-Abelian Paffien FQH state at 
v = l/2 28 ~ 32 ). This is consistant with the fact that those 
four Z2A topological states all have 3 degenerate ground 
states on torus and do not have time reversal symme- 
try. Therefore the four Z2A topological states are non- 
Abelian states with excitations that carry non-Abelian 
statistics described by Ising topological quantum field 
theory. 

We believe that the four Z2A topological states de- 
scribed by {Cfc} = 0111, 1011,1101, 1110 arc also non- 
Abelian states described by Ising topological quantum 
field theory. However, 0111, 1011,1101, 1110 states are 
different from the 1000, 0100, 0010, 0001 states since 
thay have different ground state degeneracy on (oo) lat- 
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tices. Some concrete constructions of those non-Abelian 
topological states and other topological states discussed 
in this paper are given in appendix B. 

We also like to point out that the mean-field Hamil- 
tonian is a Hamiltonian for a superconductor. The Z 2 
topological indices (k provide a classification of trans- 
lation invariant 2D topological superconductors as dis- 
cussed in Ref. 33. 



VII. CONCLUSION 

In this paper, we study topological phases that have 
the translation symmetry using the projective approach. 
We concentrated on a class topological phases described 
by Z2A PSG. We find that 2D Z2A topological phases can 
be further divided into 16 classes, which can be described 
by four Z 2 topological variables (fc = 0, 1 at the four k = 
points. Through the projected SC wave-functions, we 
obtain the topological degeneracies and the crystal mo- 
mentum for those 16 classes of topological states. This 
allows us to identify the topological phase of the Wen's 
plaquette model as the 0110 Z2A topological order. In 
addition, it is predicted that massless fermionic excita- 
tions appear at the quantum phase transition between 
different topological orders with translation invariance. 

This research is supported by NSF Grant No. DMR- 
0706078, NFSC no. 10228408, and NFSC no. 10874017. 



APPENDIX A: DEFINITION OF 
ANTI-COMMUTATION FRAMES 

In this appendix we define anti- 
commutation bases. An anti-commutation base 
(M{ a/3 }, M{ a ipty, M{ a iif)ii}, ...) is a maximum set for 
several 4x4 matrices which anti-commute each other 

{%},%«} = 0, 
{M {a , , } ,M {a „ „ } } = 0, 
{M {aP} ,M w ,p, l} } = 0, 

If a matrix can be decomposed into an anti-commutation 
base, 

U = /Z C {"/3} M {«/3} 

{«,/?} 

one has the determine of the matrix U as 

det U = [c\ afi} + c\ al[jl} + c 2 {a „p„ } + ...} 2 . 

Any matrix out of the base Mr & py can be consid- 
ered as a particular direction in the base, M, & ^y — 
J2{ a 0} C{ Q /3}-W{a/3}- When one add it to U, we have the 
determine as 

det U = [{c {a[j y + c {a[j y) 2 + {c{ a >pi} + c {a >p>}) 

+ (c{ a "p"} + £{ Q "/3"}) 2 + (Al) 



A famous example is five 7 matrices, M^y, Myi 3 y, 
M{ny, M{- 2} and M{ 01 y making up of an anti- 
commutation base. If the matrix has the form as 

u = zZ C {<*P} M {<*P}, {a/3} = {33}, {13}, {11}, {02}, {01}, 

{«,/?} 



one has the eigenvalue of U as 



-{33} '-{13} Til} 

and the determine of U as 



detU= ( c| 33} + c\ 13} + c\ 11} 



C {02} + C {01} 



"{02} "t" c {01} 



M{ 30 }, M{ 12 } and My 2 2} make up of an anti- 
commutation base, I. If there is a matrix as 



U 



C{ a p}M {al3}l 



{a/3} = {30}, {12}, {22}. 



{«,/3} 

One has the eigenvalue of U as 



±Jc 



{30} + C {12} 



+ C 



{22} 



and diagonalize U = c {3fi yM {3fi y + c {1 , 2 }M {1 ^y + 
C{2,2}^{2,2} into a 4 x 4 matrix as 



{30} ^ °{12} 



{22} 



/ 1 \ 
0-100 
10 

V -1 / 



-^{33} 7 -^{31} an d M{ 02 } make up of another anti- 
commutation base, II. One can diagonalize U = 
C{3,3} M {3,3) + C{3,i} M {3A} + C{o,2}^{o,2}, into another 
4x4 matrix as 



/ 1 \ 
0-10 
-10 

\0 1/ 



,2 

'{33} 



-{31} 



-{02} 



APPENDIX B: THE EXAMPLES OF THE 
ANSATZS OF DIFFERENT Z2A TOPOLOGICAL 
ORDERS 

In this part we give one example for each type of topo- 
logical order. The ansatzs U{k) in momentum space for 
the 16 classes topological orders are given by 



0000 : 



0101 : 



Ci 
Ci 



B v 



0011 : 



0110 : 



B x 



B x 



+y 



10 



1111 



d 
-Ci 



1100 : 



a 3 
—B x 



1000 : 



a 3 
C 3 



0010 : 



a 3 
C 4 



1010 



1110 



0111 : 



0001 : 



a 3 
-B n 



a 3 



o -c x y 



a 3 



-c. 



a 3 
Ci 



1001 : 



1011 



1101 : 



0100 : 



u 



—B 



x+y 



The parameters above are defined as 



a 3 
-Co 



a 3 
-d 



a 3 
C 2 



<r } = 

B x = 

By = 

B x +y 







1 + \ cos k x + \ cos k y u 

—I — j cos k x — j cos fc y 



cos fc x i sin /c^ , 
—i sin fcj; — cos k x J ' 

cos ky i sin fcj, \ 
—i sin fcj, — cos k y J ' 

/ cos(k x + ky) isin(k x ■ 
y —i sin(k x + k y ) — cos(fc x 



-ky) 

hky) 



Ci - 

c 2 = 

c 3 = 

c 4 = 



' cos k x — \ co%{k x 
sin k x — ? 

' cos k x — \ cos(k x 

sin k x — i sin k y 



ky ) COS ky 



fc„ Sill Kj X T t SHI Ky 

\ cos(k x + k y ) 



cos k x + 



cos(k x 

sin k x — i sin /c y 



+ fc y )+cosfc !y sink x +ii 
sin /c y — cos k x + ^ < 

— cos fcc 

/c^, "I - ^ sin ky 

1" ^y) ^ COS fcy 



COS(fc ;r i -i ; 

sin k x — i sin /c y 

1 cos fcy 



sin k x t i 
k x 2 cos{k x -+- fcy 



sin /c y 

— cos fcy 



+ \ cos fc y sin k x -\- % sin 

cos k^. cosiykrg ~\~ k^ 

L 



i sin k y 

'X ~t~ ky ) + cos k v 



sin k x + i sin A, y 
— cos fca; — cos(k x + k y ) + | cos k y 



One can check above results by the following table 

(0,0) (0,tt) (tt,0) (tt.tt) 

1 + 2 cos k x + j cos /c y 

COS fca; 

cos /c y 

COs(/Ca; + k y ) 

COS fcr + COS k y — i COs(/c.j. + fcj,) 
COS fc x — COS k y — I COS^ + fc y ) 
COS + COs(/Ca; + k y ) + | COS ky 

cos + cos(k x + ky) — 2 cos k y 



> 


> 


> 


> 


> 


> 


< 


< 


> 


< 


> 


< 


> 


< 


< 


> 


> 


> 


> 


< 


< 


> 


< 


< 


> 


< 


< 


< 


> 


> 


< 


> 
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